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We study the magnetic excitation spectra of resonant inelastic x-ray scattering (RIXS) at the 
L-edge from undoped cuprates beyond the fast collision approximation. We analyse the effect of the 
symmetry breaking ground state on the RIXS process of the Heisenberg model by using a projection 
procedure. We derive the expressions of the scattering amplitude in both one-magnon and two- 
magnon excitation channels. Each of them consists of the isotropic and anisotropic contributions. 

The latter is a new finding and attributed to the long range order of the ground state. The presence 
of anisotropic terms is supported by numerical calculations on a two-dimensional spin cluster. We 
express the RIXS spectra in the form of spin-correlation functions with the coefficients evaluated 
on the cluster, and calculate the function in a two dimensional system within the 1/5 expansion. 

Due to the anisotropic terms, the spectral intensities are considerably enhanced around momentum 
transfer q = 0 in both one-magnon and two-magnon excitation channels. This finding may be 
experimentally confirmed by examining carefully the q-dependence of the spectra. 

PACS numbers: 


I. INTRODUCTION 

Resonant inelastic x-ray scattering (RIXS) has at¬ 
tracted much interest as a useful tool to investigate ex¬ 
cited states in solids^. The L-edge RIXS experiments 
have been recently carried out with high energy resolu¬ 
tion in transition-metal compounds, which have revealed 
magnetic excitations as spectral peaks in the low-energy 
regiork^Tii. Starting from the undoped cuprates, the ac¬ 
tivity spreads over, rapidly and widely, doped high-Tc 
cuprates^^— , nickelates^, pnictides^, 5 d transition metal 
compound o^^d^ , and so on. Among them, the investi¬ 
gation on cuprates is one of the most active fields due 
to the relation with high-T,. superconductivity. Stimu¬ 
lated by these experiments, theoretical efforts to eluci¬ 
date mechanism of the magnetic RIXS in cuprates also 
have developedi^^— . But rich information involved in the 
L-edge RIXS data such as momentum and energy trans¬ 
fer dependence as well as polarization dependence never 
ceases to require further reliable and convincing theories 
more than ever. 

The L-edge resonance in undoped cuprates is described 
by the second-order dipole allowed process that a 2p-core 
electron is prompted to an empty orbital by ab¬ 

sorbing photon and then an occupied 3d electron com¬ 
bines with the core hole by emitting photon. When the 
3d orbital in the photo-emitting process is different from 
the one in the photo-absorbing process, the excitations 
within the 3d orbitals are brought about, which are called 
as d-d excitatioi*i^. When the 3d orbitals in the photo¬ 
absorbing and photo-emitting processes are the same 
x^—y^ orbital but their spins are different, magnetic exci¬ 
tations with spin flip could be generated^Si^i. Even if the 
spins are the same, the spin-conserving excitations could 


be brought about by the presence of the core hole during 
its finite lifetimei^i^— . This process could be described 
only when it is treated beyond the fast collision approxi¬ 
mation (EGA) that no relaxation could take place in the 
intermediate state because of the short Iifetime2£i2i. 

In our previous paper o^^d^ , we have analysed the pro¬ 
cess leading to the final states in the second-order pro¬ 
cess, and have clarified how the spin excitations are taken 
place around the core hole site beyond the EGA. In one 
dimension, the analysis has been straightforward, since 
the spherical symmetry in spin space remains intact in 
the ground state, while in two dimensions under the 
antiferromagnetic ordered phase, the analysis has been 
rather complicated due to the breaking of spherical sym¬ 
metry. In both cases, we have obtained the scattering 
amplitudes in an invariant form with respect to the po¬ 
larization vectors of the incident and scattered x-rays, 
and spin operators. Disregarding possible effects of the 
symmetry breaking ground state, we have obtained spin 
excitations extending to neighbours of the core hole site. 
Such excitations have been clearly observed in a one¬ 
dimensional system GaGu203^ and in two-dimensional 
systems Sr2Gu02Gl2^ and La2Gu04^. 

However, under the presence of the antiferromagnetic 
long-range order, it may be reasonable to presume that 
the scattering amplitudes include anisotropic terms as¬ 
sociated with the direction of the staggered moment, 
since the second-order process could be affected by the 
anisotropy originated from the broken symmetry of the 
ground state. This observation contrasts to neutron scat¬ 
tering, in which the scattering amplitude is directly de¬ 
scribed by the interaction Hamiltonian between the spins 
of neutron and electron. The purpose of this paper is to 
clarify the presence of anisotropic terms in the scatter- 
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ing amplitude under the presence of the spin long range 
order by analysing the second-order process on a model 
of undoped cuprates, where the low-lying excitations are 
described by the Heisenberg model. In the scattering am¬ 
plitudes summarised in an invariant form, we obtain the 
anisotropic terms, which include a vector characterizing 
the staggered moment. 

To estimate quantitative impact of the anisotropic 
terms, we evaluate them by carrying out numerical anal¬ 
ysis on spin clusters. For a cluster of 13 spins, which is 
regarded as a model of two-dimensional cuprate, various 
terms in the scattering amplitudes are calculated. We 
verify the anisotropic terms have finite contributions. If 
the connection of the anisotropic terms with the symme¬ 
try breaking of the ground state is intrinsic, the weights of 
anisotropic terms are expected to increase with increas¬ 
ing antiferromagnetic long-range order parameter. This 
anticipation is confirmed by the numerical calculation on 
a ring of 12 spins with varying the external staggered 
magnetic field, which is given in Appendix [Cl 

Collecting up such amplitudes from all the Cu sites, 
we derive the RIXS spectra represented by spin correla¬ 
tion functions. When we investigate the correlation func¬ 
tions, analysis on a larger system might be preferable. 
This is because the spin excitations propagate through 
the entire crystal in the final state. Thus, we employ the 
I/S' expansion to the spin operators^!, which practically 
enables us to treat an infinite system. As a result, we 
can express the RIXS spectra in terms of the correla¬ 
tion functions of one-magnon and two-magnon contribu¬ 
tions. Since two magnons are excited close to each other, 
their mutual interaction is important. We treat multi¬ 
ple scattering of two magnons by following the method 
previously developedi^. It turns out that the correlation 
functions for both one-magnon and two-magnon channels 
have anisotropic contributions in addition to isotropic 
ones. We hnd that the anisotropic terms produce sub¬ 
stantial enhancement on the RIXS intensities for momen¬ 
tum transfer q close to the F point in both channels. This 
shows a sharp contrast to the fact that the contributions 
from the isotropic terms vanish at q = 0 in both chan¬ 
nels. Experimentally, polarization analysis may help to 
clarify the existence of the anisotropic terms, since the 
polarization dependence is completely different between 
the one-magnon and the two-magnon spectra. 

The present paper is organized as follows. In section 
im we describe the second-order dipole allowed process 
responsible for RIXS process. In section IIIIl we anal¬ 
yse the RIXS process paying attention to the influence 
of the symmetry breaking on the scattering amplitude, 
in which anisotropic terms are derived in an invariant 
form. In section IIYI we evaluate numerically the am¬ 
plitudes of creating excitations on a finite size of two- 
dimensional cluster under a molecular field on the bound¬ 
ary. In section |Vl we derive the RIXS spectra in terms of 
the spin-correlation functions, which are treated with the 
1/5 expansion to spin operators. The RIXS spectra con¬ 
sisting of one-magnon and two-magnon excitations are 



(a) (b) (c) 


FIG. 1: A schematic illustration of the RIXS process at the 
I/ 3 -edge of cuprates. Filled and open circles stand for the 
states occupied by electrons and holes, respectively. Red and 
blue arrows denote the magnetic moments of 3d holes. Wavy 
green and pink dotted arrows represent the photon and the 
transitions, respectively, (a) An incident photon excites a core 
2p3 electron into the empty 3d state. The site of the excited 
electron is chosen as the origin of the crystal-hxed coordinate 
system with x, y, and a axes (or also called as a, b, and 
c axes). The direction of the unit vector of the staggered 
magnetic moment at the origin is called as Bm, which dehnes 
the spin coordinate system with x', y', and z' axes. The 
spin quantization axis is coincident with z' axis, (b) In the 
intermediate state, spin degree of freedom vanishes at the 
core-hole site, (c) The excited electron recombines with the 
core-hole by emitting the scattered photon. 

calculated. Section EH is devoted to the concluding re¬ 
marks. In Appendix El absorption coefficients at the 
L2— and La-edges are briefly discussed. A short com¬ 
ment on the projection procedure is given in Appendix 
IH In Appendix o we show how the expansion coeffi¬ 
cients develop with increasing the staggered moment in 
a finite-size ring under the external staggered field. Ap¬ 
pendix [D] outlines the 1/5 expansion in the Heisenberg 
model. 


II. SECOND-ORDER DIPOLE ALLOWED 
PROCESS 

We briefly explain L-edge RIXS in cuprates (See hgure 
[T]). The RIXS process at the copper L-edge may be de¬ 
scribed by the electric dipole (El) transition between the 
2 p-core states and the 3 d states. The 2 p states are split 
into two levels with the total angular momentum j = ^ 
and which are discriminated as L2 and L3 edges, re¬ 
spectively, due to the strong spin-orbit interaction. Be¬ 
cause each Cu atom has one hole in the orbital 

in undoped cuprates such as La2Cu04 and Sr2CuCl202, 
we employ a hole picture. Then, the Ei transition may 
be expressed by the interaction between photon and hole 
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as, 


Hint = W 


lint =W ^ 

q,a,2,m,cT 


D°‘{jm,a) 

\/2a^ 




difjQ 


H.C., 


( 2 . 1 ) 

where Cqa annihilates a photon with four-vector q = 
(q, Wq) and polarization a. The hi represents the cre¬ 
ation operator of the 2p hole with jm at site i, and dia 
denotes the annihilation operator of the 3d hole with the 
orbital and spin a at site i. The w is a constant 
proportional to r^i? 3 d(r)i? 2 p(?’)dr where Rsdir) and 
i? 2 p(c) are the radial wave-functions for the 3d and 2p 
states of Cu atom. The D°‘{jm,a) describes the depen¬ 
dence of the El transition amplitude on the 2p core-hole 
angular momentum and the spin of the 3d hole. 

In the El transition at the L 2 - and La-edges, the ini¬ 
tial photon having q = qi, a = ai excites 2p core hole 
into empty 3d state, which decays back into the 2p state 
by emitting the final photon having q = qt^ a = a{. 
The RIXS spectra associated with this process may be 
expressed as 


interaction J only between the nearest neighbour sites. 
In the thermodynamic limit, the ground state of idmag on 
a square lattice is spontaneous symmetry broken phase, 
that is, long-range ordering antiferromagnetic phase. 

We write the ground state \g) of idmag as 

\9) = \ml) + \mt), ( 3 . 2 ) 

where | f) and | |) represent the spin states at the origin, 
and IV'o) and IV'o) are constructed by the bases of the rest 
of spins. We assume that a core hole is created at the 
origin as a result of absorbing photon (figure [T|). In the 
intermediate state, the spin degrees of freedom is lost at 
the core-hole site, since the 3d hole in the x^—y^ orbital is 
annihilated by the 2p-3d dipole transition. Note that the 
Hamiltonian in the intermediate state is similar to that 
for a system with a non-magnetic impurity introduced 
into antiferromagne t^^i^^ . 

Denoting as the eigenstate of the intermediate 
Hamiltonian with eigenvalue we can express the 
second-order amplitude in ( 12 . 21 ) as 


W{q{a{-,qiai) = 27r^ 

/' 

X S{Eg -|- Wi — Efi — uif), (2.2) 

with qi = (qi,Wi), qi = (qf,a;f), |$i) = 4 iailg)|vac), 
1 $/') = |/')|vac), where |g) and |/') represent the 

ground state and excited states of the matter with en¬ 
ergy Eg and Efi, respectively. Note that f refers to the 
final state of the photon and /' refers to the excited state 
of the electron. The |vac) is the vacuum state for pho¬ 
tons. The eigenstate and its energy of the intermediate 
state are referred to as |n) and En, respectively. 

Incidentally, since the final state in the absorption co¬ 
efficient A{uj) is the intermediate state in the RIXS, we 
have 

A{uj)=2TrY,\{n\Hintm\^S{Eg+L0-En). (2.3) 

n 

The explicit form is summarised in Appendix lAl 


'f’lHir 


rj \ /.n I TX. 


III. MAGNETIC EXCITATIONS AROUND THE 
CORE-HOLE SITE 

In undoped cuprates, the low-energy excitations may 
be well described by the two-dimensional antiferromag¬ 
netic Heisenberg Hamiltonian on a square lattice, 

ifmag = T ^ S, . S,q (3.1) 

where Si denotes the spin one half operator at site i, and 
{iA') indicates that the summation runs over nearest- 
neighbour pairs. Since our focus is not on a discussion of 
the magnetic dispersion, we have adopted the exchange 


Hint\n){n\Hint\g) na,/- t ■ l\ 

2_^ , I F _ F — ^ AJm,a) D Ajm,a) 


g 




with 


R{e') = 


uj{-\- — ec 


.. + iT-€'J 


(3.4) 


where Cg represents the ground state energy of Hmag. 
The Ccore denotes the energy required to create the 3d^°- 
configuration and a 2p core hole in the state \jm}. The T 
stands for the life-time broadening width of the core hole; 
r ^ 0.3 eV at the Cu L 3 edge. Notice that the scattering 
amplitude ( 6 ) and those investigated in the remaining 
of this section are originated from the excitation of the 
single electron at the origin. A whole scattering intensity 
will be given by collecting up the amplitudes from all Cu 
sites. 

In the scattering amplitudes leading to those excited 
states, we seek the invariant form with the polarization 
vectors a.i and ccf of the incident and scattered x-rays, 
spin operators Si, and the unit vector of the staggered 
moment em. To this end, it is convenient to consider 
a general situation that Om is pointing to an arbitrary 
direction, which is denoted as axis z'. Then, for spin op¬ 
erators of the id electron, coordinate frame of x', y', z' 
axes is prepared (See figure [1] (a)). On the other hand, 
the 3d orbitals as well as spin and orbital of the core hole 
are described in the crystal-fixed coordinate frame with 
X, y, and z axes. Since the definition of the spin coordi¬ 
nate system and that of the crystal-fixed system are in¬ 
dependent, we can relate them by any method which can 
describe the transformation from the one to the other. 
We adopt here the rotation of the Euler angles a, /3, and 








4 


7 as the transformation from the crystal-fixed to the spin 
coordinate system^Z. Our final formulae do not depend 
on the specific choice of the Euler angles. 

The D^{jm^ a) in this definition is given in table I ofi^. 
Then we introduce ttf,ai) and Pa^\j]ai,a\) by 

af, ai),(3.5) 

m 

af, ai),(3.6) 

m 

where —a represent I and f for a =t and respectively. 
The Pa^\j;a{,ai) and Pi^^(j; af, ai) correspond to the 
spin-conserving and the spin-flipping processes, respec¬ 
tively. Polarizations of x-ray are along the x, y, and z 
axes defined in the original crystal axes. Since the follow¬ 
ing analysis is confined to the Pa-edge, we fix j = | and 

omit the argument in the expressions of af, ai) 

and Pi^^ {j ; af, tti). We list all the non-zero values of them 
for j = § below. 


pW(x,x) = pW(j/,y) = A, (3.7) 

pW(a:,y) = -pW(j/,x) =-sgn(a)^ cos/3, (3.8) 

15 

Pa^Hx,y) = -Pi^)(2/,x) = ^exp[i7Sgn(cr)]sin/3, 

(3.9) 

where sgn(CT) gives -1-1 and — 1 for cr =t and I, respec¬ 
tively. Note that P(°)(af,ai) and Pa^\a{,ai) are zero if 
a\ = z and/or ai = z. This results from the fact that the 
process is restricted with the hole of the — y^ orbital 
in the ground state. 


A. Scattering channel with changing polarization 


the core-hole site, we approximate |P) by a linear combi¬ 
nation of the states = S^lg) and = A=*=|g), 

where X = | Sj with j running over the nearest 
neighbour sites around the core-hole site. The number 
of the nearest neighbour sites z is four and two for two 
and one dimensions, respectively. Spin raising and lower¬ 
ing operators on the core-hole site and neighbouring site 
are defined as = Sq ± iS^ and ± iX^ , 

respectively. 

Since the |'0i^^)’s are not orthogonal to each other nor 
normalized, we need to introduce the density matrices 
to project \F) onto A 

procedure to determine the expansion coefficients is given 
in Appendix [B] where the projection formalism is utilised. 
It may seem strange the non-orthonormal set is used in 
the expansion. However, since the procedure described 
in Appendix |B] can determine the expansion coefficients 
uniquely for the finite number of the projected states, the 
non-orthonormal set can have a one-to-one correspon¬ 
dence with some orthonormal set, for instance, by means 
of Gram-Schmidt process. Since the physical meaning of 
each element of the non-orthonormal set is much clearer 
than that of the orthonormal one, we use the former. 

Then, the final state is approximately expressed as 

|F) « ^P(i)(af,a0 ^ 

<7 


This expression is rearranged as 


\F) ^ Pl^>{ai,aO /o\7(c.0^o“lg)+ClMX 


(1). 


f(i)/ 




f(i) 


(1)/ 


(3.12) 


The coefficients are given by 


As seen from (EH) and (13.91) . both Pi°^(af,ai) and 
Pi^^(af,ai) have off-diagonal elements with ai and a{. 
This implies that the scattering channel with changing 
photon polarization includes both the spin-flipping and 
spin-conserving processes. Let us investigate them sepa¬ 
rately in the following. 


1. Spin-flipping process 

The final state arising from the spin-flipping process 
may be written as 

( 3 - 10 ) 

a r} 

Assuming the magnetic excitation associated with the 
creation of core-hole at site 0 has a local character around 


1 = 1,2 

X X I (3-13) 

r) 

1 = 1,2 

X (i/XlXlt)l</>.)^«)(</>.IV'o^)- (3-14) 

r) 

The can be constructed from by replac¬ 
ing (/3^^^“^)i,j with ■ Let us examine each 

coefficients appeared in (13.121) . We suppose that the core 
hole site belongs to ‘up spin’ sublattice. This does not 
mean 5'()’|(?) = 0 when \g) is the symmetry broken an¬ 
tiferromagnetic ground state. That is, the spin can be 
raised even at the ‘up spin site’. Then, for example, 
if spin-flip excitation takes place at the core-hole site. 








5 


two channels, from up spin to down spin and vice versa 
should be survived. Each cahnnel experiences different 
surroundings in the intermediate state through the sec¬ 
ond order process, which is materialized due to the fact 
that the core-hole has a finite life-time. As a result, both 
channels aquire different values of the coefficients. Sim¬ 
ilar explanation is also valid for the spin-flip process at 
the nearest neighbour sites. 

In the presence of the antiferromagnetic long-range or¬ 
der, the coefficients /q and for a =t are 

expected to be different from those for a =4-. Then, let 
us divide them into two parts as follows. 

= /o^^(wi)+sgn(CT)A^^^’o(wi), (3.15) 

-bsgn(cr)A|L\\(wi). (3.16) 


utilises the density matrix pij = {'ipil'ijjj). Hence the final 
state in this channel is approximately expressed as 

a- ij 

X XI (3-19) 


This relation is rewritten as 


\F') 


X afu) 


5o^^(wi)So|| 




(wi)X|| 


|g), 

(3.20) 


It has been confirmed that Aj_.o(wi) = Aj_^i(a;i) = 0 and 
in ttie absence of the long-range order 

for one-dimensional systen>i^. Therefore, the A^^^(a;i) 
stands for the anisotropic part of the coefficient, while 
fn^\uji) represents the isotropic part of coefficient. Note 
that the anisotropic coefficient Aj__o(wi) and the coeffi¬ 
cients for the excitations on neighbouring sites /}^j(a;i) 
would not come out in the FCA, since the relaxation 
process in the intermediate state is disregarded. Insert¬ 
ing (13.151) and (13.161) into (13.121) with the help of (13.9L 
we notice that (13.121) with the Euler angles a, /3, 7 con¬ 
stitute an invariant form (see (3.18) for isotropic 
terms). The result is given by 


\F) - ^(Cii X af^) • {/('^ccOSox + 

Xx.o(^i)em X Soj. -I- AX(wi)em X X_L ||g), 

(3.17) 


where ciij, and Q;f_L, respectively, are polarization vec¬ 
tors of the incident and scattered photon, which are pro¬ 
jected onto the a-b plane. Operators Soj_ and Xj_, respec¬ 
tively, are Sq and X, which are projected onto the plane 
perpendicular to the direction of the staggered magnetic 
moment. 


2. Spin-conserving process 

According to (13.31) . the spin-conserving process may be 
written as 


where So|| and Xy, respectively, represent the vector oper¬ 
ators of So and X parallel to the direction of the staggered 
magnetic moment. Note that the amplitude associated 
with \'ipi) is omitted. The definition of the expansion co¬ 
efficient gn \u!i) is inferred from the projection procedure 
in Appendix |B1 We have already confirmed that ( 7 Q^^(wi) 
and ^^^^(wi) were equivalent to /^^^(wj) and re¬ 

spectively, in the absence of long-range order— There¬ 
fore, it is natural, in the presence of long-range order, to 
write them as 

= /o^^(wi) - iA[|^’(a;i), (3.21) 

- iX||^i(a;i). (3.22) 


Here A|j^Q(a;i) and A||^j(a;i) correspond to the anisotropic 
contributions of the coefficients. 

Combining the spin-conserving term (13.201) to the spin- 
flipping term (13.171) . we finally have 


|F)i = \F) + \F’) 

~ ^ “fJ-) ■ - iAX(Wi)em X Sq 


IAm Q(cUi)eni(eni * So) 


2 i 


+ X afj.) • - iA)L^_\(a;i)em x X 


( 1 ) 
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“i^iu(‘^i)em(em • X) |g). 


(3.23) 


(3.18) 

(7 r} 

where the off-diagonal elements with the polarizations 
are used for (af, ai). We approximate |F') by a linear 
combination of the states \tpi) = |g), \'ijj 2 ) = Sq |g), and 
\ 1 lj 3 ) = |g). Since these states are not orthogonal to 

each other nor normalized, we repeat the analysis that 


The terms containing e„i represent the effect of the long 
range order, that is, that of the broken symmetry in spin 
space. If em is defined on the A sublattice and the same 
em is used on the B sublattice, A||^p(wi) and A l^(wi), 
respectively, take the same value in both sublattices. On 
the other hand, the values of A||^Q(a;i) and A||^\(a;i) in 
sublattice B, respectively, are obtained by changing en¬ 
tire sign of those in sublattice A, respectively. 






















B. Scattering channel without changing 
polarization 
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Since only has the non-zero diagonal ele¬ 

ments with Oi and Of, dSSl) may be expressed as 

\F), = ^pW(a,a)|a)^ |</.,)P(e;)(0,l4-). (3.24) 

(7 rj 

We see that the FCA could not give rise to spin ex¬ 
citations in this process because the diagonal element 
(a, a) is independent of tr. Since the total spin is 
conserved, |F )2 niay be expressed by \g), Sq |g), |g), 

Sq |g), and -|-S'^X+)|g). Similar to the pro¬ 

cedure resorted in the preceding subsection, |P )2 is ap¬ 
proximated by a linear combination of these states with 
the help of the density matrix. Hence |P )2 is approxi¬ 
mately expressed as 

\F)2 - ^(ai^.afi)/f)(a;i)So-X|g) 

+ ■ Q;fi)[A(^)(wi)(em ■ So)(em ■ X) 

+ X|j p(wi)(em • So) -I- A||^j(wi)(eni • X)]|g()3.25) 

where the amplitude associated with |g) is omitted. The 
terms containing e„i represent the effect of broken sym¬ 
metry in spin space. The expansion coefficients for Sq |g) 
and X^'lg) are denoted as A||^Q(u;i) and A|^(wi), respec¬ 
tively, while those defined for Sq X^'\^ and -|- 

5'()'A+)|g) are divided into the isotropic term 
and anisotropic term A(^)(a;i). If em is defined on the A 
sublattice and the same em is used on the B sublattice, 
A(^)(wi) take the same value in both sublattices. On the 
other hand, the values of A||^Q(a;i) and A||^j(a;i) in sub¬ 
lattice B, respectively, are obtained by changing entire 
sign of those in sublattice A, respectively. 



FIG. 2: A cluster of 13 spins used to evaluate the coefficients. 
The spin at site 0 is annihilated in the intermediate state. 
Spins on the boundary are subjected to the molecular held 
from spins outside the cluster. 


in the ground state and in the intermediate state. Ta¬ 
ble U shows the calculated results at oJi giving the maxi¬ 
mum absorption coefficient with T/J = 2.4 and 1.0. The 
values for T/J = 2.4 may correspond to La 2 Cu 04 and 
Sr2Cu02Cl2. Note that the coefficients have dimensions 
of (energy)”^ as seen from right-hand side of (13.31) . The 
coefficients not shown there are small, and will be ne¬ 
glected in the calculation of the RIXS spectra in the next 
section. 

As seen from (13.231) and (13.251) . it is obvious qualita¬ 
tively that the origin of the anisotropic terms, which in¬ 
clude the unit vector representing the staggered moment 
(em), is attributed to the broken symmetry of the ground 
state in spin space. In quantitative sense, the magnitudes 
of such terms are expected to develop as increasing the 
staggered moment. This is confirmed in Appendix ICl for 
a finite-size ring of spins. 


IV. EVALUATION OF THE COEFFICIENTS 


Various coefficients defined in the preceding section 
could be evaluated by diagonalizing the Heisenberg 
Hamiltonian on finite-size clusters. Since the excitations 
are localized around the core-hole site, the calculation 
on small clusters may give reliable estimates to the coef¬ 
ficients. We consider a cluster of 13 spins shown in figure 
[2 A complication is that analysis on finite-size cluster 
cannot provide spontaneous symmetry breaking ground 
state. In order to break the spherical symmetry in spin 
space, we assume that the spins on the boundary are sub¬ 
jected to the molecular field, — J|(S'q )|, per bond. The 
expectation value of Sq is determined self-consistently 
as {Sq ) = 0.394. The coefficients in the RIXS process 
are evaluated by diagonalizing the Hamiltonian matrices 


V. ANALYSIS OF RIXS SPECTRA FROM 
UNDOPED CUPRATES 

Now, we are in a position to calculate the RIXS spec¬ 
tra. It is preferable to treating a larger system since the 
spin excitations propagate through the entire crystal in 
the final state. Thus, we employ the results of the 1/S 
expansion to the spin operators, which practically corre¬ 
sponds to taking into account of an infinite system effect 
as well as the interaction among the magnetic excita¬ 
tions. In doing so, we proceed the analysis by dividing 
the RIXS spectra into two channels, with and without 
changing photon polarization. 
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TABLE I: Various coefficients in units of 1/J in the two-dimensional cluster of 13 spins: (a) coefficients for isotropic terms and 
(b) coefficients for anisotropic terms. The incident photon energy cui is set to give the maximum absorption coefficient. 

(a) Isotropic coefficients (b) Anisotropic coefficients 

T/J 

2.4 (-0.087,-0.367) (0.177,-0.320) (-0.093,0.048) (0.071,-0.002) 

1.0 (-0.263,-0.656) (0.082,-0.891) (-0.209,0.334) (0.182,-0.207) 


A. Scattering channel with changing polarization 

Since ctfj, and ai_L are polarization vectors projected 
onto the a-b plane, ctfj, x OLi± is parallel to the c axis. 
In undoped cuprates such as La 2 Cu 04 and Sr2Cu02Cl2, 
the staggered magnetization aligns along the (1,1, 0) di¬ 
rection in the Cu02 plane^^. Therefore the anisotropic 
terms proportional to em could not come out. We col¬ 
lect up the remaining amplitudes from all Cu sites, where 
(13.231) is multiplied by the weight exp(iq- ri) at the core¬ 
hole site Yi with momentum transfer q = qj —qf. Thereby 
we obtain 

W{qi,a{;qi,oci) 

= (ai-L X af_L)^r(^)(wi;q,a;).(5.1) 

with w = Wi — Wf is defined by 

(5.2) 

where 

Z«(cui;q) = /W(a;0K'(-q) + 5,“'(-q)] 

+ - *5'^ (-q)]- (5-3) 

Here the time dependent operator of an arbitrary opera¬ 
tor A is defined as A{t) = The Fourier 

transforms of the spin operators are given by 

So(-q) = (2/A^)5^Siexp(iq-ri), (5.4) 

zGA 

Sb(-q) = (2/iV)5^Siexp(iq-ri), (5.5) 

ieB 

where the sum is taken over site i on the A or B sublat¬ 
tices. The x', y', and z' axes are defined as directing to 
(0,0,1), (1,—1,0), and (1,1,0), respectively. The spin- 
flip excitations on the neighbouring sites to the core hole 
are neglected, because their amplitudes are quite small. 

We expand the spin operators by means of magnon 
operators in the 1 /^'-expansion method, which is briefly 
summarised in Appendix [D1 In their expressions, mo¬ 
menta are defined within the first magnetic Brillouin 
zone (MBZ). When momentum q lies outside the first 
MBZ, Sa(—q) and Sf,(—q) are replaced by Sa([—q]) and 
sgn( 7 G)Sf,([—q]), respectively, where q is put back into 
the first MBZ by a reciprocal lattice vector G. That is. 


q = [q] -I- G with [q] lying inside the first MBZ. The 
sgn( 7 k) denotes the sign of 7 k, where 7 k = ^(cosfc^ -I- 
cos ky) with k in units of 1/a (a is the lattice constant). 
For example, 70 = —1 for G = (tt, tt). With these no¬ 
tations together with the magnon operators ck[_q] and 

/3[-q]: ^(^^(wqq) is expressed as 

Z^^\uji-q) = VTSMiwqq) + sgn(7G)/3j_qj) H-, 

(5.6) 

where 

M{ujqq) = ^|/(|^^(wi)[l-sgn(7G)x[q]] 

+sgn(7G)a:[q]]|. (5.7) 

The definitions of and Xq are found in Appendix [Pl 
[ (ID10|) ] and use has been made of the relations £_q = £q 
and a:_q = Xq. Therefore F(^)(cui; q, w) consists of the 
(5-function peak, which is located at 

CJ = JSZ ^1 + £[q], (5.8) 

where A = 0.1579 is the first order correction in the 1/S- 
expansion (see Appendix IdI)^. 

Figure [3] shows y(^)(wi;q, w) as a function of w for q 
along the symmetry directions with F/J = 2.4. A no¬ 
table aspect is that the intensities diverge at q = (0, 0) 
and (7r,7r). The corresponding integrated intensity is 
given by 

/(i)(a;i;q) = J y(^)(a;i;q,a;)^ = 2(2S')|M(a;i;q)|^ 

(5.9) 

Figure[3]shows q) for q along symmetry directions 

with T/J = 2.4. The intensities are enhanced around q = 
(0,0) as l/|q|. Since the contribution from the isotropic 
term vanishes around q ~ 0, the enhancement is due to 
the finite value of the anisotropic coefficient A|j^Q(a;i). 
On the other hand, the divergence around q = (tt, tt) 
is brought about by the isotropic term, which is why 
the behaviour is irrelevant of presence of the anisotropic 
term. 

It has been observed in the RIXS experiments^Ti^ that 
the intensity of magnon peak increases significantly with 
q —>■ 0. Such increase is consistent with the effects of the 
anisotropic terms. So far, the increase of intensity has 













first MBZ. See figure 7 ini^ for such RIXS spectra. 
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FIG. 3: q, oi) as a function of uj for q along symmetry 

directions. The uJi is set to give rise to the peak in the absorp¬ 
tion spectra. J = 130 meV and F/J = 2.4. The intensities 
diverge at q = 0 and (tt, tt). 


B. Scattering channel without changing 
polarization 

In order to calculate the RIXS intensity in this channel, 
we collect up the amplitudes from all the Cu sites with 
the use of (|3.25|) . We obtain 

W[qf,OLf]qi,a.{) 

= (ctii • af_L)^r^^)(a;i;q,w).(5.10) 

The correlation function q, w) is defined by 

(5.11) 

where 



FIG. 4: Intensity 7^^^(cui;q) of the d-function peak arising 
from the one-magnon excitation as a function of q along sym¬ 
metry directions. The uJi is the photon energy giving rise to 
the peak in the absorption spectra. F/J = 2.4. The (black) 
solid and (red) broken lines are the results with and without 
including A^p(a;i), respectively. 


been interpreted simply as the contribution from elastic 
scattering. To confirm the effects of anisotropic terms, 
it may be necessary to examine carefully the spectra 
with subtracting systematically the contribution of elas¬ 
tic scattering around q ^ 0. 

It should be noted here that there exists non-linear 
terms which make the one-magnon excitation split into 
three-magnon excitations in the second order correc¬ 
tion of the 1 /S'-expansioni2i^. Accordingly Y^^'> (wi; q, w) 
contains the energy continuum of the three-magnon ex¬ 
citations in addition to the 5-function peak mentioned 
above. The contribution from the three-magnon excita¬ 
tion grows gradually when q is near the boundary of the 


Z(2)(a;i;q) = /f (a;i)[(S • X)„(-q) + (S • X)b(-q)] 

+ A(2)(a;0[(.5^'X^')J-q) + (5^'X^')b(-q)] 

+ ••• . (5.12) 

The Fourier transform {S^ )o(—q) and 

( 5 'M )&(—q) for = x', y', and z' are introduced as 
follows. 

(5^'X'^')J-q) = {2/N)i J E 

ieA 5 

= (2/iV)i^5,^'(k)5r'(-[k + q])7[k+q], (5.13) 

k 

(5^'x^'),(-q) = i2/N)i Y. E 

zgB i 5 

= (k)5^'(-[k + q])7[k+q]Sgn(7G). 

k 

(5.14) 

Here the sum over 5 is carried out on the nearest neigh¬ 
bour sites around site i. 

Expanding Z^^')(uji;q) in terms of magnon operators 
within the l/5'-expansion(see Appendix iDl). we obtain 

Z(^'>(uJi;q) = (2S)Y^(^i'’^>^)(^-[q+k]/^l + ' " > (5-15) 


with k running within the first MBZ, and 


A^(a;i;q, k) 

= {-(1 -b 7q) [a;k + sgn(7G)a:[q+k]] 


-b [7k + sgn(7G)7[q+k]] [1 + sgn(7G)a;[q+k]Xk]} 


- a(^H ^0 


^[q-kk ]4 


(1 + 7q) [ik + sgn(7G)a;[q+k]](5.16) 
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FIG. 5: q, tj) as a function of cu for q along symme¬ 

try directions. The cJi is set to give rise to the peak in the 
absorption spectra. J = 130 meV and F/J = 2.4. 

This expression is valid even when q is outside of the first 
MBZ. Note that when A(^^(a;i) = 0, N(uJi;q,k) vanishes 
at q = (0, 0) and (tt, tt)^. Note also that the isotropic 
terms of the two-magnon part are the same as those ob¬ 
tained for the Ff-edge RIXS, where no anisotropic term 
exists22r^. 

From (I5.15L we see that q, uj) consists of the 

energy continuum of the two-magnon excitations. Since 
two magnons are created at neighbouring sites through 
x-ray scattering, inclusion of the magnon-magnon inter¬ 
action is crucial to obtain the spectral shape. As al¬ 
ready discussed in^^, the magnon-magnon interaction in 
the 1 /S'-expansion could be divided into a separable form 
so that the t-matrix of the scattering is neatly evalu¬ 
ated. We resort to the similar evaluation. Figure [5] shows 
(wi; q, w) as a function of uj for q along the symmetry 
directions. We find rapid enhancement of the intensity is 
brought about by the presence of the anisotropic terms 
as |q| goes to (0,0). Without them, in contrast, the in¬ 
tensity diminishes in this limit as shown in figure 8 ofi^. 
We see the peak energy decreases with |q| approaching 
zero. At q = 0, the peak energy becomes very close to 
zero, ^ 0.025eV. It may be a difficult task to distinguish 
the spectral peak from the elastic peak. However a care¬ 
ful study on the q-dependence of the spectra may clarify 
such effect of the anisotropic terms. 

The frequency integrated intensity /*'^^(a;i;q) may be 
given by 


FIG. 6: Frequency-integrated intensity 7^^1(a;i;q) of the con¬ 
tinuous spectra arising from two-magnon excitations as a 
function of q along symmetry directions. The cJi is the pho¬ 
ton energy giving rise to the peak in the absorption spectra. 
F/J = 2.4. The (black) solid and (red) broken lines represent 
the results obtained with and without the anisotropic term, 
respectively. 



FIG. 7: Schematic view of the scattering geometry. The an¬ 
gle between the incident and scattered x-rays is fixed at 130 
degrees. The scattering plane contains the c axis. 


at a region where q is away from (0,0). 


C. Polarization dependence 


/(2)(a.i;q) ^ J (a.q q, c)^ 

= (2^)2^|iV(a;i;q,k)^ (5.17) 

k 

Figure |6] shows /^^^(wi;q) for q along symmetry direc¬ 
tions. Notice that {uji; q) diverges logarithmically 
when |q| approaches zero. It demonstrates a possibil¬ 
ity that this logarithmic enhancement can be recognized 


We consider a scattering geometry used in the exper¬ 
iments of La 2 Cu 04 - and Sr2Cu02Cl2^. It is schemati¬ 
cally shown in figure [7] for q along (0,0) — (0, tt) direction, 
where the angle between the incident and the scattered x- 
ray is fixed at 130 degrees. The scattering plane includes 
the b{x) and c(z) axes. 

Then ccij^ = (1,0,0) for the a polarization and ccij, = 
( 0 ,0^,0) for the tt polarization in the incident x-ray, 
while Q!f_L = (1,0,0) for the a' polarization and a.fx = 
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(0,af^,0) for the tt' polarization in the scattered x-ray. 
Thereby the RIXS spectra may be expressed as 

W{qf,af,qi,oti) 

( (aj^)2r(i)(a;i;q,a;), a ^ n 

_ 2 \ I r(2)(a;i;q,a;), a ^ a' 

4wiWf V15/ I 7r-)>cr' 

[ {aJaf)'^Y^'^\u!i;q,uj), tt-)> tt' 

(5.18) 

The one-magnon term and the two-magnon term 
are separated by the polarization. Accordingly the 
polarization analysis is useful to clarify the contribution 
of Y^^\ For other directions of q, we could obtain the 
similar formulas separated by polarizations. 


VI. CONCLUDING REMARKS 


Appendix A: Absorption coefficient 


Since the 2p core hole is quite localized in real space, 
the absorption coefficient is well approximated by the 
sum of the intensities on each lattice site. Therefore, af¬ 
ter averaging the polarization, the absorption coefficient 
A(j, Wi) at the Cu L 2 - and La-edges may be expressed as 


A(j, wi) oc - ^ I—-- 

^ [wi + eg - e. 




_e/12_pr2’ 

v„j -r -L 


(Al) 


where Ccore depends on j . By substituting the eigenvalues 
and the eigenstates evaluated on finite-size clusters into 
(El), we obtain A(j,Wj). For F comparable or larger 
than the excitation energy e'^, A{j,uji) is close to the 
Lorentzian curve. The calculated A(j, Wj) for F/J = 2.4 
in a two-dimensional Heisenberg model has been shown 
in figure 3 ini^. 


We have studied the magnetic excitations on the L- 
edge RIXS from undoped cuprates beyond the FCA. Em¬ 
phasis is on how the symmetry breaking of the ground 
state affects the magnetic RIXS spectra. It is found that 
the spin excitations are brought about at neighbouring 
sites in addition to the core-hole site. We have shown 
that the anisotropic terms emerged in the scattering am¬ 
plitudes as a direct consequence of the broken symmetry. 
The fact contrasts sharply with the case in neutron scat¬ 
tering, where the amplitude is described through the in¬ 
teraction Hamiltonian between the spins of neutron and 
electron. The presence of such anisotropic terms has been 
supported by the calculation on a one-dimensional finite- 
size ring of spins under the staggered external field and 
on a two-dimensional cluster with the molecular field act¬ 
ing on the boundary. Collecting up such amplitudes on 
all the Cu sites, we have expressed the RIXS spectra in 
the form of spin correlation functions, which have been 
calculated within the I/S'-expansion. The anisotropic 
terms have made the RIXS intensity considerably en¬ 
hanced as |q| goes to zero. Such enhancement could be 
confirmed experimentally by observing carefully the spec¬ 
tra around q = 0. With a little further improvement on 
energy resolution, it would make the distinction possible 
as achieved in Sr 2 lr 04 , in which band-splitting, predicted 
by a theory^, conspicuous around q = (0, 0) had been 
discerned by recent experiment^. We believe our present 
emphasize on the anisotropic terms originated from the 
antiferromagnetic long range order might be insightful 
when one analyses the systems with short range order 
such as the doped high-Tc cupratesi^^— . 
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Appendix B: Projection onto non-orthogonal bases 

We try to project a state \F) onto the non-orthogonal 
states |4>i)’s as 

= (Bi) 

i 

Operating (^'ij from the left side of (IBII) . we obtain 
ft = where = (4'j|F) and pij = (4'i|4'j). 

Therefore the expansion coefficients /i’s are given by 
fi = As long as the number of the pro¬ 

jected states remains finite, this procedure uniquely de¬ 
termines the expansion coefficients. 


Appendix C: One-dimensional ring of spins under 
the staggered external field 

We examine how the anisotropic terms develop in con¬ 
cert with the development of the staggered moment. This 
is achieved easily in one-dimensional system, since there 
is no long range order in the ground state. To control 
the staggered moment, we apply the staggered external 
field. Thereby the Hamiltonian of the system is given by 

ifmag = J ^ S, • S, + Llex (Cl) 

dj) * 

where the field strength is denoted as i/ex- We consider a 
system made of 12 spins of S' = ^ with periodic boundary 
condition, as shown in figure [51 

Since the total-spin component along the z' axis is zero 
in the ground state, the Hamiltonian matrix is repre¬ 
sented by a matrix with 924 x 924 dimensions. Diagonal¬ 
izing the Hamiltonian matrix, we obtain the ground-state 
wavefunction. The inset in figure [5] (a) shows the stag¬ 
gered magnetization (Sg ) as a function of i/ex, which 
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FIG. 8: A ring of 12 spins (S' = i) used to evaluate the 
coefficients for the spin excitations in RIXS. The spin at site 
0 is annihilated in the intermediate state. 


TABLE II: Coefficients for isotropic terms in units of 1/J in 
the one-dimensional ring of 12 spins (S = |). The incident 
photon energy uji is set to give the maximum absorption co¬ 
efficient. 

r/J /f^(a;i) 

2.0 (-0.013,-0.500) (0.043,-0.025) (0.366,-0.219) 

1.5 (0.025,-0.672) (0.088,-0.051) (0.633,-0.365) 

1.0 (0.050,-1.009) (0.154,-0.125) (1.101,-0.894) 


increases with increasing iLex- The intermediate state, 
on the other hand, consists of 11 spins, since the spin 
degree of freedom is lost at the core-hole site. Thereby 
the Hamiltonian matrix in the intermediate state is rep¬ 
resented by a matrix with 462 x 462 dimensions in the 
subspace of Using the eigenvalues and 

eigenfunctions for the Hamiltonian of the intermediate 
state together with the ground state, we evaluate (I3.11L 
(13.1911 . and others for the coefficients of the spin excita¬ 
tions. 

By setting Hex = 0, we first evafuate the isotropic 
terms in the absence of the anisotropic terms. Table m 
shows the coefficients for isotropic terms for several val¬ 
ues of r/ J with oji being fixed at the value to give the 
maximum absorption coefficient. The values T/ J = 2.0 
and 1.5 may correspond to CaCu203^ and Sr2Cu03^, 
respectively. The coefficient for the spin-fiip 

excitation on neighbouring sites is much smaller than 
/o^^(wi). The coefficient f 2 ^\uji) for the So ■ X term is 
comparable to the coefficient /g^^(wi) for the spin-flip 
term. It grows with decreasing T/ J, as was discussed 
iiiH. 

Next, we turn our attention to the anisotropic terms. 
FiguresEKb) and (c) show the absolute values of the coef¬ 
ficients as a function of staggered moment for T/J = 2.0. 
They demonstrate that the anisotropic terms develop 
with increasing staggered moment as expected. Note 
that the magnitudes of the isotropic terms vary gradually 
and slightly diminish rather than increase with increasing 



FIG. 9: Various coefficients measured in units of J~^ as a 
function of the staggered moment evaluated in a ring of 12 
spins. T/J = 2.0, and uui is set to give rise to the maxi¬ 
mum absorption coefficient, (a) The (black) solid and (red) 
broken curves correspond to |/q^^ | and 1 / 2 ^^ |, respectively, (b) 
The (black) solid and (red) broken curves represent and 

|A^^1|, respectively, (c) The (black) solid, (red) broken, and 
(blue) dotted curves show |A||^g|, lAy^jl, and |A||^g|, respec¬ 
tively. The inset in Panel (a) shows the staggered moment as 
a function of Hex/ J- 


staggered moment as shown in figured (a). 


Appendix D: l/S'-expansion 


Here, we briefly summarise an introduction of the 1/5 
expansion. The emphasis is on the definitions of the 
quantities used in the main text. The details are rele¬ 
gated to the references such asi^. Assuming two sublat¬ 
tices in the antiferromagnetic ground state, we express 
spin operators by boson operators as^l 



= 5-aJai, 

(Dl) 


= (5-)t = v/^/,(5)a„ 

(D2) 


= -S + b]hj, 

(D3) 

St 

= (5-)t = V^blf.iS), 

(D4) 
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where and bj are boson annihilation operators, and Here, 




(D5) 


with ni represents a]at and for £ = i and j, re¬ 
spectively. Indices i and j refer to sites on the A and B 
sublattices, respectively. Using (jDll) - (ID4l) . is ex¬ 

panded in powers of I/S', 

= -JS^Nzl2 + + • • • , (D6) 

where JV and z are the number of lattice sites and that 

(n) 

of nearest neighbour sites, respectively, ffmag stands for 
the n-th order term in the 1/S-expansion. The Fourier 
transforms of the boson operators are introduced within 
the first MBZ, 

a, = (2/N)iJ^ ake''^"^ , b, = {2/N) ^ ^ 6ke‘''"^ . (D7) 

k k 


Then, with the help of a Bogoliubov transformation, 

al = £kal + mk,5_k, 6-k = rukal -I- 4/3-k, (D8) 

we could diagonalize iSmig as 

= JSzJ^(ek-l)+JSzJ^ek(alak+/3l/3k). (D9) 

k k 


]/ 2ek^ “ a.k4,(D10) 

= a/i -7k. 7k = 

^ S 

where S connects the origin with the nearest neighbour 
sites. The expression for iFmag becomes^^ 

^mlg = ek(a];ak + 

_ 7C~ 

+ <5G(ki-fk2-k3-k4) 

ki.k2.k3,k4 

X 4i 424344 (4a^^/3lj.^/3_k2<ak3£3kik2k3k4 "I ) . 

(D12) 

with A = ^ X]k(l ~ Ck)^- For the square lattice, 
A = 0.1579. The Kronecker delta (5 g( kid-k2 —kg —k4) in¬ 
dicates the conservation of momenta within a reciprocal 
lattice vector G. In the second term of (IDI2I) . only the 
relevant term representing scattering of two magnons is 
shown explicitly. The vertex function in a symmet¬ 
ric parametrization as well as omitted terms are found 
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